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Abstract 

We report measurements of the branching fractions for B° tt+tt", 
K+TT-, K+K- and ir°7r°, and B+ 7r+7r°, K+n^, K'^tt+ and K+W. The 
results are based on 10.4 fb"-*^ of data collected on the T(4S) resonance at the 
KEKB e+e" storage ring with the Belle detector, equipped with a high mo- 
mentum particle identification system for clear separation of charged tt and 
K mesons. We find B{B^ tt+tt') = (0.56 Igjg ± 0.04) x 10-^ ^(5° ^ 
i^+TT-) = (1.93 t'oM t^ol) X 10-^ B{B+ K+ttO) = (1.63 t^H) x 
10-5, B{B+ ^ K^Tr+) = (1.37 toll tofs) x 10~^ and ^(5° ^ K'^tt'^) = 
(1.60 Iggg Iq 2?) ^ , wheiB thc first and second errors are statistical and 
systematic. We also set upper limits of B(B+ — tt+tt'^) < 1.34 x 10"^, 
B{B^ K+K-) < 0.27 X 10-^ and B{B+ K+K°) < 0.50 x 10"^ at the 
90% confidence level. 

PACS numbers: 13.25.Hw, 14.40.Nd 



3 



The charmless hadronic B decays B vrvr, Ktt and KK provide a rich sample to test 
the standard model and to probe new physics Of particular interest are indirect and 
direct CP violation in the tctt and Kir modes, which are related to the angles 02 and 03 of 
the unitarity triangle, respectively Measurements of branching fractions of these decay 
modes are an important first step toward these CP violation studies. However, experimental 
information is rather limited, and the only published results come from one experiment . 
One of the key experimental issues is the particle identification (PID) for separation of the 
high momentum charged vr and K mesons. This is one of the primary reasons that the B 
factory experiments have been equipped with specialized high momentum PID devices. 

In this paper, we report the first results of the Belle experiment on charmless hadronic 
two-body B decays into vrvr. Kit and KK final states. The decay modes studied are tt'^tt~, 
K+TT-, K+K- and K^tt^ for 5° decays, and n+n^, K+n°, K°Tr+, K+W, for B+ decays. 
For the modes with K^ mesons, only K^ tt+tt^ decays are used. Throughout this paper, 
the inclusion of charge conjugate states is implied. The results are based on data taken by 
the Belle detector at the KEKB asymmetric e+e^ storage ring [Q. The Belle detector 
is equipped with aerogel Cerenkov counters (ACC) configured for high momentum PID. 
The data set consists of 10.4 fb~^ data taken at the T(4S) resonance, corresponding to 11.1 
million BB events, and 0.6 fb~^ data taken at an energy ~60 MeV below the resonance, for 
systematic studies of the continuum qq background. 

Primary charged tracks are required to satisfy track quality cuts based on their impact 
parameters relative to the interaction point (IP). Kg mesons are reconstructed using pairs 
of charged tracks that have an invariant mass within ±30 MeV/c^ of the known K'^ mass 
and a well reconstructed vertex that is displaced from the IP. Candidate vr" mesons are 
reconstructed using 7 pairs with an invariant mass within ±16 MeV/c^ of the nominal 7r° 
mass. The B meson candidates are reconstructed using the beam constrained mass, rribc = 
\/-^beam ' Pbi ^ud the energy difference, AE = Eb - -Ebeam, where -Ebeam = v^/2 ~ 5.290 
GeV, and ps and Eb are the momentum and energy of the reconstructed B in the T(4S) 
rest frame, respectively. The signal region for each variable is defined as ±3cr from its central 
value. The resolution in rribc is dominated by the beam energy spread and is typically 2.7 
MeV/c^. The AE resolution ranges from 20 to 25 MeV, depending on the momentum and 
energy resolutions for each particle. Normally we compute AE assuming a vr mass for each 
charged particle. This shifts AE downward by 44 MeV for each charged K meson, giving 
kinematic separation between the /itt"*" and hK~^ [h = tt, K) final states. In modes with 
7r° mesons, both the m^c and AE distributions are asymmetric due to 7 interactions in the 
material in front of the calorimeter and energy leakage out of the calorimeter. We accept 
events in the region rribc > 5.2 GeV/c^ and |Aii^| < 0.25 GeV for the h~^h~ and K^h^ modes, 
and —0.45 < AE < 0.15 GeV for the h^n^ and Kgir^ modes. In this kinematic window, the 
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area outside the signal region is defined as a sideband. The signal reconstruction efficiencies 
after the kinematic window cut are 65% for h^h~ , 33% for Kgh~^, 50% for /i+vr", and 24% for 
K^Ti^, according to a GEANT [§ based Monte Carlo (MC) simulation. The MC tracking 
efficiency is verified by detailed studies using high momentum tracks from D, rj and K* 
decays. The reconstruction efficiencies for high momentum Kg and vr" mesons are tested by 
comparing the ratio of the yield of K^tt^ to K^ii^n^ and —>■ K~tt~^7[^ to 

— > K~7i~^, respectively, between data and MC simulation. From these studies, we assign 
a relative systematic error in these efficiencies of 2.3% per charged track, 12% per Kg and 
8.5% per 7r° meson. 

The background from b ^ c transitions is negligible. The dominant background is from 
the continuum qq process. We suppress this background using the event topology, which is 
spherical for BB events and jet-like for qq events in the T(4S) rest frame. This difference 
can be quantified by using several variables including the event sphericity, S, the angle 
between the B candidate thrust axis and the thrust axis of the rest of the event, 6^, and 
the Fox- Wolfram moments ifj = J^ij IPillPjlPiicosOij), where the indices i and j run over 
all final state particles, pi and pj are the momentum vectors of particles i and j, Pi is the 
l-th Legendre polynomial, and % is the angle between particles i and j. We can also use 
the B fiight direction, 9b, and the decay axis direction, 9hh, which distinguish BB from qq 
processes based on initial state angular momentum. 

We increase the suppression power of the normalized Fox- Wolfram moments, Ri = 
Hi/ Ho, by decomposing them into three terms: Ri = R^'+Rf°+R^° = {H^' +H^°+H°°) / Hq, 
where the indices ss, so, and oo indicate respectively that both, one, or neither of the parti- 
cles comes from a B candidate. These are combined into a six term Fisher discriminant 
called the Super Fox- Wolfram |TD[ defined as SFW = Ef=i (a«^f° + PiRD, where ai and 
/3; are Fisher coefficients and /=2,4 for ai and The terms i?f* and RfZi 3 are excluded 
because they are strongly correlated with m^c and AE. In the h~^h~ modes, for example, 
SFW gives a 20% increase in the expected significance compared to R2. 

We combine different qq suppression variables into a single likehhood, £^(^5) = Y[i'^l{gg), 
where the -^^^(95) denotes the signal (gg) likelihood of the suppression variable i, and select 
candidate events by cutting on the likelihood ratio TZg = Cs/{C-s + Cqq). For h^h" and 
Kgh'^, the likelihood contains SFW, cos 9b, and cos9hh- In modes with 7r° mesons, the qq 
background is significantly larger. In this case, we first make a loose cut on cos9t- Next, we 
extend SFW to include cos 9t and S, and form the likelihood using this extended SFW and 
cos 9b- In each case, the signal probability density functions (PDFs) are determined using 
MC simulation and the qq PDFs are taken from m^c sideband data. The performance of 
TZs varies among the modes with efficiencies ranging from 40% to 51% while removing more 
than 95% of the qq background. The tt'^h^ mode calls for a tighter cut with an efficiency 
of 26%. The error in these efficiencies is determined by applying the same procedure to the 
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D^TT^, —>■ K~iT~^ event sample and comparing the cut efficiencies between data 
and MC. The relative systematic error is determined to be 4%. 

The high momentum charged vr and K mesons (1.5 < ph± < 4.5 GeV/c in the laboratory 
frame) are distinguished by cutting on the '7i{K) likelihood ratio 71tt(k) = '^■K{K)/{'^n + ^^.k), 
where Ctt{k) denotes the product of each tt^K) likelihood of their energy loss (dE/dx) in 
the central drift chamber and their Cerenkov light yield in the ACC. Each likelihood is 
calculated from a PDF determined using MC simulation. The PID efficiency and fake rate 
are measured using vr and K tracks in the same kinematic range as signal, with kinematically 
selected D*~^ D^tc^, K~tt^ decays. The efficiency and fake rate for vr mesons are 

measured to be 92% and 4% (true vr fakes K), whereas those for K mesons are 85% and 
10% (true K fakes vr), respectively. The relative systematic error in the PID efficiency is 
2.5% per charged it or K meson. 

Figure || shows the m^c and AE distributions in the signal region of the other variable, for 
the n^TT^ , K^TT^ and Kgir'^ modes. Each m^c and AE distribution is ffited to a Gaussian 
signal plus a background function. The rribc and AE peak positions and rubc width are 
calibrated using the — * D^n^, K^n^ data sample. The AE Gaussian width 

is calibrated using high momentum K^tt^ and —>■ K^tt^ decays. The mbc 

background shape is modeled by the ARGUS background function |Tl| with parameters 
determined using positive AE sideband data. A linear function is used to model the shape 
of the AE background; the slope is fixed at the value determined from the m^c sideband. The 
signal yields are determined from the AE fits where there is kinematic separation between 
the /iTT"*" and hK'^ decays. The 7t~^tt~ and K^7t~ fits include a component to account for 
misidentified backgrounds. The normalizations of these components are free parameters. 
The extracted yields are listed in Table |. The cross-talk among different signal modes is 
consistent with expectations based on PID fake rates. No excess is observed in the K^K~ 
and K^Kg modes. 

Figure § shows the rribc and AE projections for the tt+tt", K^h^ and K^tt'^ modes. For 
these modes, since the AE distribution has a long tail, a two-dimensional fit is applied 
to the iribc and AE distributions. The signal distribution is modeled by a smoothed two- 
dimensional MC histogram, while the background distribution is taken to be the product of 
the and AE background functions discussed above. The signal and background shapes 
are determined following the same procedure as for the h^h~ and K^h^ modes. The AE 
resolution is calibrated using K^tt^ti^ decays where the 7r° is reconstructed in the 

same kinematic range as the signal. For the tt+tt" mode, since the cross-talk from K'^n'^ is 
expected to be large and the AE separation is less than la, the K~^7i^ component is fixed 
at its expected level. The obtained yields are listed in Table |I[ 

The systematic error in the signal yield is determined by varying the parameters of the 
ffiting functions within ±la of their nominal values. The changes in the signal yield from 
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each variation are added in quadrature. These errors range from 1% to 6%. In the K^tt~ 
mode, the AE background normahzation is influenced by an excess around —175 MeV. In 
this region, we expect to observe a few background events from B decays such as -B — pvr, 
K*TT, and i^*7 (for modes with 7r° mesons), based on a MC simulation [0,|l3 in all signal 
modes. To estimate their effect, we either exclude the negative AE sideband from the fit 
or add these components to the fit based on MC histograms. The resulting change in the 
signal yield, ranging from 4% to 10%, is added in quadrature to the above systematic error. 

Table | summarizes all results. The statistical significance (S) is defined as 
\J —'^ ln(£(0)/£max), where £max and £(0) denote the maximum likelihood with the nominal 
signal yield and with the signal yield fixed at zero, respectively [|12|. The final systematic 
error is the quadratic sum of the relative error in the signal yield (iVs), the reconstruction, 
PID, and continuum suppression efficiencies, and the number of BB pairs (1%). If E < 3, 
we set a 90% confidence level upper limit on the signal yield {N'^'^') from the relation 
/o " C{Ns)dNs/ C{Ns)dNs = 0.9, where C{Ns) denotes the maximum likelihood with 
the signal yield fixed at A^^- The branching fraction upper limit (U.L.) is then calculated by 
increasing N^'^' and reducing the efficiency by their systematic errors. 

In summary, using 11.1 million BB events recorded in the Belle detector, the charge 
averaged branching fractions for B tt+tt^, K^tt^ , K^tt^, K^tt^, and K^tt'^ are measured 
with statistically significant signals. For the '7r"'"7r° mode, an excess is seen with marginal 
significance. No excess is observed for the K^K^ and K^K^ modes. For these modes, 
90% confidence level upper limits are set. The results are listed in Table |. In Table 0, 
we list some ratios of branching fractions based on these measurements. Recent theoretical 
work [ij suggests that the ratio B{B^ 71^71^) /B{B^ tt+tt") is relevant for extracting 
02, the ratio B{B^ K^tt^)/B{B^ K^n^) is relevant for determining the contribution 
from electro- weak penguins, and the remaining four ratios are useful to constrain ^3. All 
the branching fraction and ratio results are consistent with other measurements Our 
results confirm that B{B^ K^n^) is larger than B{B^ tt+tt^), and indicate that 
B{B^ h^Ti^) and B{B^ K^tt^) seem to be larger than expected in relation to the 
B^ h^n^ and B^ K^tc^ modes based on isospin or penguin dominance arguments [|l|. 
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(GeV/c^) AE (GeV) 



FIG. 1. The nibc (left) and AE (right) distributions, in the signal region of the other variable, 
for B ^ a) vr+vr^, b) K^tt^ and c) Kgir^. The fit function and its signal component are shown by 
the solid and dashed curve, respectively. In the 7r"'"7r~ and K~^Tr~ fits, the cross-talk components 
are shown by dotted curves. 




(GeV/c^) AE (GeV) 



FIG. 2. The rubc (left) and AE (right) projections for B ^ a) 7r+7r°, b) K+tt^ and c) K^tt^. 
For if+7r°, a K mass is assumed for the charged particle. The projection of the two-dimensional fit 
onto each variable and its signal component are shown by the solid and dashed curve, respectively. 
In the 7r+7r° fit, the cross-talk from if+7r° is indicated by a dotted curve. 
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TABLE I. Summary of the results. The obtained signal yield (Ns), statistical significance (S), 
efficiency (e), charge averaged branching fraction (B) and its 90% confidence level upper limit 
(U.L.) are shown. In the calculation of B, the production rates of B^B~ and B^B^ pairs are 
assumed to be equal. In the modes with mesons, Ng and e are quoted for K^, while B and 
U.L. are for . Submode branching fractions for — > 7r"'"7r~ and 7r° — ^ 77 are included in e. 
The first and second errors in Ng and B are statistical and systematic errors, respectively. 

Mode Ns S 6 [%] B [xlO'^] U.L. [xlO-^] 



50^ 




177 +7.1 +0.3 

J^'-' -6.4 -1.1 


3.1 


28.1 


0.56 tofo ± 0.04 




B+- 




-,r, 4 +5.1 +1.2 
-4.3 -1.6 


2.7 


12.0 


n 70 +0.38 +0.08 
-0.32 -0.12 


1.34 


B^^ 




f^r, q +10.6 +2.7 


7.8 


28.0 


1 QO +0.34 +0.15 
^^•^•3 -0.32 -0.06 




B+ ^ 




nA Q +7.6 +0.6 
'^^■•^ -7.0 -2.0 


7.2 


19.2 


1 f-n +0.35 +0.16 
-•^•""J -0.33 -0.18 




B+ ^ 




in +4-3 +0.4 
^^•'^ -3.6 -0.1 


3.5 


13.5 


1 07 +0.57 +0.19 
-"-•■J' -0.48 -0.18 






► i^V 


A +3.8 +0.4 
-3.1 -0.6 


3.9 


9.4 


1 fin +0-72 +0.25 
^■^^ -0.59 -0.27 






K+K- 


A 9 +3.8 
^•^ -0.2 




24.0 




0.27 


B+^ 


K+W 


^■^ -0.0 




12.1 




0.50 



TABLE II. Ratio of charge averaged branching fractions {B) for B tttt, and Kir decays. The 
first error is statistical and the second is systematic. The correlation and cancellation of systematic 
errors are taken into account. A 90% confidence level upper limit in the first ratio, is calculated 
using a similar method as the upper limit of B described in text. 



Modes Ratio 

B{B+ tt+ttO) /B{B^ TT+TT-) < 2.67 

2B{B+ ^ if+7rO)/^(SO ^ K+tt') 1.69 f'oll 

B{B^ 7r+7r-)/B{B^ ^ K+n') 0.29 t'oH 

B{B° ^ i^+7r-)/2 BiB^ ^ K^tt^) 0.60 1°:^ t^H 

2 B{B+ ^ K+Tr^)/B{B+ ^ K\+) 2.38 tUl t^H 

B{B' K+n-)/BiB+ K'^-^) 1.41 t^H 
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